The strongly interacting Quark Gluon Plasma, 
and the critical behaviour of QCD at imaginary 
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We explore the highly non-perturbative hot region of the QCD phase diagram close to Tc by 
use of an imaginary chemical potential fj, which avoids the sign problem. The number density 
and the quark number susceptibility are consistent with a critical behaviour associated with the 
transition line in the negative fj? half-plane. We compare the analytic continuation of these results 
with various phenomenological models, none of which provides a satisfactory description of data, a 
failure on which we make some comments. These results complement and extend the information 
obtained via the analysis of the susceptibilities evaluated at zero ^, yielding a simple description of 
the candidate strongly interacting QGP phase. As a byproduct of our analysis we investigate the 
Polyakov loop and its hermitian conjugate. Our data offer a vivid evidence of the importance of the 
complex nature of the functional integral measure, which results in L{ii) 7^ L{y?j for a real chemical 
potential ^. 

PACS numbers: 12.38 Gc, 11. 15. Ha, 12.38.Mh 



I. INTRODUCTION 



Theoretical arguments and experimental evidence sug- 
gest that hadronic matter undergoes a transition to a 
plasma of quarks and gluons at high temperature 
At extremely high temperatures quarks and gluons are 
nearly free, and should be described by the Stefan- 
Boltzmann law with the appropriate degrees of freedom. 
When temperature is not much larger than the critical 
temperature - say, Tc < T <~ 2Tc - strong interactions 
among the constituents give rise to non-perturbative ef- 
fects. In short, at large T the QGP is a gas of nearly 
free quarks, which becomes strongly interacting at lower 
temperatures T = (1 - 3) 0, Q. 

Several proposals have been made to characterise the 
properties of the system in such non-perturbative phase. 
For instance the above mentioned strong interactions 
might be enough to preserve bound states above T^, while 
coloured states might appear^ deeply affecting the ther- 
modynamics of the system l4|. Analytic techniques are 
being refined more and more, so to be able to capture the 
features of the system closer and closer to Tc [5] . Model 
theories of quasi particle physics have been considered as 
well [6, 7]. In this work we study this interesting dynam- 
ical region by lattice QCD simulations at T ~ 1.1 Tc (the 
reason for this choice will be clear in the following), and 
a nonzero baryon density. 

In principle, lattice QCD simulations at non-zero 
baryon density are plagued by the sign problem [8| . How- 
ever, it has been realised that this problem can be cir- 
cumvented thanks to physical fluctuations, which grow 
relatively large in the Quark Gluon Plasma phase. In 
this work we adogted the ima gina ry chemical potential 
approach [1, ^ [ll|, [H, [H [11 111 , which avoids the 
sign problem and makes it possible conventional Lattice 
QCD simulations. The interested reader mi ght want to 
consult refs. [13, [3 for recent reviews and [l9l [20| for 
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FIG. 1: Schematic phase diagram for four flavor QCD in the 
T, plane. The candidate sQGP phase is bound by the chiral 
(pseudo) critical line in the negative half-plane 

more pedagogical introductions into the subject. 

Other studies in the quark gluon plasma phase have ad- 
dressed the higher temperature regime 21, 22, 23, 24, 25|. 
Here we analyze in detail the non-perturbative behaviour 
close to Tc, in the candidate strongly coupled Quark 
Gluon Plasma (sQGP) region (some preliminary results 
have appeared in [26]). We note that in the sQGP re- 
gion the chiral critical line lies in the imaginary chemical 
potential plane, and that such a chiral line ends in the 
proximity of the endpoint of the Roberge- Weiss line [131 . 
We focus our analysis on the particle number and its sus- 
ceptibility, on the chiral condensate, and on the Polyakov 
loop, and we find that the results are consistent with 
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those expected of a critical behaviour associated with 
the critical line at imaginary chemical potential. Hence, 
the numerical results are compatible with simple power 
law behaviour of the equation of state as a function of 
the imaginary chemical potential fij , yielding a modified 
form of the Stefan-Boltzmann law. 

The rest of this paper is organised as follows. Section 
II is devoted to the analysis of the particle number and 
the chiral condensate, which are related by the Maxwell 
equation. In Section III we discuss the behaviour of the 
Polyakov loop. It turns out that our approach offers 
a particularly simple description of an apparent puzzle, 
and, at the same time, gives a direct evidence of the phase 
of the determinant at nonzero, real chemical potential. 
The implication on the equation of state are summarised 
in Section III, while Section IV discusses our results in 
the light of phenomenological proposals, and alternative 
lattice approaches. Last Section is a short summing up. 
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FIG. 2; Rpifii) = n{fii)/n{fii) function of /i/ , show- 

ing a very clear evidence of a deviation from a free field be- 
haviour 



II. THERMODYNAMICS OF THE HOT PHASE 
CLOSE TO Tc 

Let us remind ourselves of the critical lines in the phase 
diagram in the T, fj,^ plane (Figure [1]): at high tempera- 
ture there is the Roberge Weiss transition at /i = ttT/S, 
associated with the phase of the Polyakov Loop, ending 
at T ~ Tnw At lower temperatures the QGP region is 
limited by a chiral transition at negative /i^, which con- 
tinues into the physical chiral transition at positive /i^, 
i.e. real chemical potential pU. [T^. 

While II approaches 7rT/3 at a constant temperature 
T ~ Tjiw the chiral transition approaches the Roberg 
Weiss transition. Within the current numerical accu- 
racy the endpoint of the two transitions cannot be re- 
solved, and the nature of the critical behaviour around 
T — Tftn^^ii — ttT/S is an interesting question in it- 
self. If T is slightly larger than T^w we are approaching 
the Roberge Weiss transition, if slightly lower we hit the 
chiral transition, and at T = Tbw we might expect inter- 
esting critical phenomena whose universality class is not 
known a priori. Note that in the chiral limit the /i = 
transition should be of first order, and we do not expect 
any tricritical point along the critical line at a real chem- 
ical potential. A possible occurrence of an endpoint at 
finite mass in the T, /i plane depends on dynamical de- 
tails, which are not known, and are not relevant for the 
present study. 

We have then carried out simulations on a 16'^ x 4 
lattice and four flavor of staggered fermions at (3 — 5.1, 
which, according to our previous results, yields T ~ T^w, 
endpoint of the RW transition,. Fermions are fully de- 
generate, with a bare dimensionless mass (a being the 
lattice spacing) to = ma = 0.05. For our lattice the 
value of the (dimensionless) ft = ^la which is relevant for 
the Roberg Weiss transition reads p, — tt/12 (the tem- 
perature being T = l/{aNt) and in our case Nt = 4). 
With a slight abuse we will omit in the following the 



hat-notation, nevertheless measuring fi and T in unit of 
inverse lattice spacing. 

First, we check our data for the particle number against 
a simple free field behaviour. We have numerically com- 
puted the free field results for real chemical potential on 
a 16"^ X 4 lattice, and niq — 0.05, and we have fitted them 
to an expression motivated by one dimensional QCD [2^ , 
which turns out to be an excellent parametrisation: 



n{li)fr 



3sinh(/VT) 
K + cosh{fi/T)' 



(1) 



yielding the free field results for the number density as a 
function of imaginary chemical potential 



n{lii)fr 



3sin(^/T) 
K + cos{n/T)' 



(2) 



We then considered the ratio between the numer- 
ical results and such free field results Rf{ij.i) = 
n{ij,j)/n{iii) free (Figure [21). We observe a clear depen- 
dence of RpifJ-i) on fij: the results are qualitatively dif- 
ferent from a free field and the discrepancy cannot be 
accounted for by any simple renormalisation of the de- 
grees of freedom. This behaviour should be contrasted 
with that of Fig. 12 of Ref. where the results at high 
temperature did differ from a free field behaviour by a 
constant factor very close to one. 

As a second attempt at interpreting our data in terms 
of the simplest parameterisations, we consider an analogy 
with the Hadron Resonance Gas model [13, HH . We will 
discuss more fully these aspects in the last Section. 

In a nutshell, the HRG describes the system as a gas of 
weakly interacting resonances. The pressure of the HRG 
model reads: 

P(T,,)-P(r,^ . ^^(T)(cosh(^)^l) (3) 



rp4 



T 



F(T) ^ I dmp{m)C^fK,C^) (4) 
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and the argument of the hyperbolic cosine, NcHq, tells us 
that in the hadronic phase one can only excite baryonic 
degrees of freedom. 

Let us remind ourselves here that general arguments 
guarantee that the partition function is periodic at imagi- 
nary /i, and that the strong coupling analysis shows that 
periodicity is smooth at low temperature. Hence, the 
number density reads 



boF sm{kNcNtfJ.i) 



(5) 



When HRG holds true, 
should suffice. 



one term in the Fourier series 



= K{T) siniN.Nt^ii) 



(6) 



A cross check with the Hadron Resonance which uses the 
Taylor approach requires the computation of an infinite 
number of derivatives [3l| , while the Fourier analysis - 
possible with the imaginary chemical potential approach 
- needs only one parameter fit. 
In Figure [3] we display the ratios 



nipi) 



shi(3fii /T) 
3sin(^//r) 



(7) 
(8) 



RsifJ-i) should be a constant for a simple hadron gas 
(cfr. again Ref. p^). while, mutatis mutandis, Rq{fii) 
should be a constant for a "hadron gas" made of quarks. 
Both Rb and Rq stay constant only for a short inter- 
val of chemical potential, indicating the region where 
n{fii) oc fj,j. The deviation from a linear behaviour (for 
/i/ > 0.05) , as well as from the simplest trigonometric 
parameterisations (for /i/ > 0.15), are evident from the 
plot. 

We now move on and propose to describe the particle 
number in the critical region as 



n{jii) = Afii{fij - fijy 



(9) 



This ansatz for n{^i) takes into account that n{^) should 
be an odd function of /i. Moreover, it reproduces a sin- 
gular behaviour of the quark number susceptibility at a 
genuine critical point. Namely, the most divergent part 
of the quark number susceptibility Xq behaves as 



1 



(10) 



where 7 = 1 — a, while /ij = 7r/12 if the critical point 
coincides with the Roberge- Weiss line. 

We then fit our data to Eq. ^ with /ij either open 
or constrained. A fit to our entire interval with un- 
constrained ^i'J gives A = -0.94(4), = 0.0804(2), 
a = 0.28(2) with a reduced = 2.4. We checked the 
stability of these results by choosing different ranges in 
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FIG. 3: Rx{l-ii) (defined in eqs. 2 and 3) as a function of fii , 
showing the limitations of a linear approximations, as well as 
the deviations from the simplest trigonometric parameterisa- 
tions motivated by a Hadron Resonance Gas model 



chemical potential, and we obtained the exponent a rang- 
ing from 0.34(8) and 0.26(3), fxf between 0.078(4) and 
0.091(12), with reduced ranging between 1.8 and 5. 
(see Figure 2] for two representative fits) . 

If we constraint /ij^ = (7r/12)^ the quality of the fits 
decreases giving a reduced ~ 12. If we limit the fitting 
interval to /i/ > 0.15, we need to add a constant to the 
function to approximate the regular component in this 
interval to obtain a reduced = 3.5, with a = 0.12(1). 
All in all, constraining fj.'j does not improve the results. 

We now go back to the quark number susceptibility 
as entailed in Eq. pUj) . for which our results indicate 
7 = 1 — q; ~ 0.7. Figure [5] shows Xg obtained by numer- 
ical differentiation of n{ni). The numerical quality is of 
course poor, but, anyway, a fit to the form of Eq. (fTU]) 
with an open /ij gives 7 = 0.66(16), while a fit with con- 
strained /ij gives 7 = 0.44(22), in agreement with the 
above estimates, within the large errors. 

The chiral condensate < ipif) > can be inferred from the 
chiral equation of state in either phases, and also in the 
presence on an explicit chiral symmetry breaking term. 
To make a closer contact with thermodynamics, we con- 
sider the Maxwell relation (the temperature is constant 
and its dependence is omitted) [3^ : 



dn{fi,m) d < ip^ > {fi,m) 



dm 



(11) 



Considering the m dependence in the expression for 

,c2/ \ ,,2\a 



n{fi, m) = F{'m)fj,{fi'j (to) — fi Y 



(12) 



where F{m) depends only on the quark mass, and using 
Eq. (fTTj) we arrive at 



< VV'> (m) - H{m){,iY{m)-ii') 



2\a+l 



(13) 



+ K{m){^xf{m) - m')"+ < ^V- >o 
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FIG. 4: n(/i) fitted to the form predicted by tlie simple critical 
beliaviour at imaginary /i Eq. ([9} 
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FIG. 6: Numerical results for the chiral condensate, and 
the results of the fit to a functional form inferred from the 
Maxwell equation Eq. 
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FIG. 5: The quark number susceptibility obtained by numer- 
ical differentiation of the results for the quark number 
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FIG. 7: The Polyakov loop fitted to the form predicted by a 
simple critical behaviour at imaginary /i Eq. (|15p . 



where H{m),K{m) depends only on the mass and 
< ipip >o is an integration constant which can be fixed 
e.g. by the chiral condensate at zero chemical potential. 

We have then fitted the chiral condensate to the lead- 
ing term at ~ 

<i>il; >= K{b- ^i"^)"" + A (14) 

obtaining a nice fit with a reduced — 0.79, A = 
0.552(6), b = 0.06628(8), K -0.63(2) and = 0.47(2) 
in reasonable agreement with a estimated from the num- 
ber density. By constraining the fitting interval /i/ > 0.2 
the sub-leading contributions are less important, and 
ttc — 0.32(12) gets even closer to a. 

It might be interesting to compare this critical be- 
haviour with that of the endpoint of QCD from model 
field theories [ssj . 



III. THE POLYAKOV LOOP 

In the same spirit we have fitted the traced Polyakov 
loop L = TrP to a power law form 

Li^lI) ^ {^lf - ^,'Jf (i5) 

Figure [7] displays the results of the fit of the absolute 
value of the Polyakov loop, which looks indeed satisfac- 
tory. 

It is interesting to look in more detail into the be- 
haviour of L = TrP. 

The Polyakov loop P satisfies the same relation as the 
quark propagator at nonzero chemical potential [38l | 

Pi^l) = pt(_^) (16) 
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This relation implies that, while both L = TrP and 
L — TrP^ are real at real chemical potential, L ^ L, a,s 
noted in m, mm mill]. We will show that the results 
at imaginary chemical potential offer a particularly sim- 
ple illustration of these ideas, as well as a direct evidence 
of the complex phase of the determinant. 

The asymmetry at real chemical potential is easily un- 
derstood by considering the distribution of the Polyakov 
loops in the complex plane: for /x = 0, since the Z3 centre 
symmetry is broken by the dynamical quarks, the root 
corresponding to the phase </> = is preferred, and the 
average is non zero. A non zero, positive chemical poten- 
tial encourages forward propagation: the distribution of 
the phases is further peaked a.t 4> — 0, while the two other 
phases have the same probability. Hence, the Polyakov 
loop remains real, and the final average is again real, 
different from zero, and slightly larger than the one at 
zero density. L instead describes backward propagation: 
again the Polyakov loop remains real, however its length 
is reduced, hence L{^) ^ L{fi). 

Notice that at a first nai've look it may sounds strange 
that, while configuration by configuration the Polyakov 
loop and its hermitian conjugate are always the complex 
conjugate of each other, their expectation values, even 
being real, differ from each other. However it should be 
clear that the complex nature of the functional integral 
measure plays an essential role in this respect, since the 
real part of the expectation value is not just the expec- 
tation value of the real part. In that sense the fact that 
L{fi) ^ L{p) is directly linked to the complex nature of 
the fermion determinant, thus also giving a qualitative 
feeling about the severeness of the sign problem. 

An apparent puzzle then arises when one considers 
the behaviour at imaginary chemical potential: there the 
measure is real and one can show that the absolute value 
of L and L are equal as well as their real parts. What is 
then the fate of the asymmetry which is present at real 
chemical potential? 

Consider 



n(ni) ' — &- 
Im Polyakov ^--a- 



(17) 



where Eq. (|16p has been used in last equality. Lo/e are 
respectively even and odd in /i. Remember that the ana- 
lytic continuation to imaginary chemical potential of an 
even function is real, while the analytic continuation of 
an odd function is purely imaginary. Hence, the analytic 
continuation of the even observable L^ip) = L(ii) -\- L{p) 
at imaginary chemical potential is the real part of L{^j)\ 
while the analytic continuation of Lo{fJ.) = L{p) — L{p) 
is the imaginary part of L{fj,i) at imaginary fi. L itself 
has no definite ^-parity and its analytic continuation de- 
velops an imaginary part. 

We conclude that we must search for the analytic con- 
tinuation of the asymmetry L ^ L which is present at real 
chemical potential in the imaginary part of the Polyakov 
loop, which is non-zero in presence of an imaginary chem- 
ical potential. Figure [5] shows the imaginary part of 
L(iii): it is different from zero, offering a clean, direct 
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FIG. 8: The imaginary part of the Polyakov loop, divided by 
the coefficient of the linear term of a fit at small /i, as a func- 
tion of the imaginary chemical potential, demonstrating the 
relevance of the phase of the determinant for a real chemical 
potential; in the same plot we show the number density, again 
normalised by the coefficient of the liner term. 



evidence of the asymmetry L ^ L, hence of the complex 
phase of the determinant at real chemical potential. In 
the same Figure |8] we have also plotted n(/i/) (both L 
and n with an appropriate normalisation), to show its 
correlation with L{fii). This correlation is in agreement 
with the lattice interpretation of the number operator 
n as 'counting' the links winding forwards minus those 
winding backwards, and it should become exact in the 
heavy quark limit when the model reduces to a Polyakov 
loop model [li,!!^. 

Let us now come back to the critical behaviour of \L\: 
since Im(i(/i/) oc n(fj,i), we might expect that |L(/i/)| 
approaches zero with a similar power law. The fit of 
Figure [7| gives an exponent of 0.45(2), in the same range 
as 7. It would be interesting, of course, the study of the 
string tension and other correlators in the same range of 
chemical potentials. 



IV. ANALYTIC CONTINUATION TO REAL 
CHEMICAL POTENTIAL 

Finally, we can analytically continue^ our results for 
n(/i/), obtaining 



i(A*)=A^(^f +m')^ 



(18) 



^ The singularity in the complex fi plane might well limit the radius 
of convergence of the Taylor expansion. However, the proposed 
- Pade' like'- parametrisation does not suffer from this prob- 
lem [id l4lj and it should not come as a surprise that we are 
able to analytically continue beyond the radius of convergence of 
the Taylor series: the analytic continuation is then valid for all 
the real values of the chemical potential fi. 



Analytic Continuation 
Stephan Boltzmann 




FIG. 9: n{fi) from analytic continuation, together with a free 
field behaviour is shown for comparison. The fits suggest that 
the slower increase observed in the interacting case with re- 
spect to the free case can be described by on overall exponent 
smaller than one 



with a ~ .3. The results are shown in Figure [9l In the 
same diagram we also plot the free field results, as an 
indicative comparison. Note that coefficient A has a non 
trivial dimension, again indicating that the system is not 
free. 

It is then amusing to notice that by using simple argu- 
ments from the theory of critical phenomena we arrive at 
a modified (lattice) Stefan-Boltzmann law, which would 
correspond to a = 1, and a large /xj^ ~ 0.5. In this 
framework a large /ij can be interpreted as a spinodal 
point at imaginary chemical far away from the Roberge 
Weiss line. 

Obviously, Eq. accounts for a slower increase of 
the particle density closer to Tc than in the free case. 
This is expected on physical grounds, as well as from 
the behaviour of the susceptibilities [H, [i^ li?!. lisLliol. 

, and of course accounts for the behaviour observed in 
Figure [H 

From the results above, we conclude that the data in 
the candidate region for a strongly coupled QCD are ac- 
counted for by a conventional critical behaviour: clearly, 
a free field behaviour would have been incompatible with 
it. In other words, the nonperturbative features of the 
plasma are closely related with the occurrence of the crit- 
ical line at negative fi^. 



V. THERMODYNAMICS AND PARTICLE 
CONTENT 

Early in this note we have contrasted our data with 
some very simple parameterisations motivated by the 
HRG approach. Here we discuss this point in some more 
detail, in particular we wish to examine the proposal that 
in the deconfined, strongly interacting region considered 



here one might observe either coloured and colourless |42l | 
bound states: in short, at large T the QGP is a gas of 
nearly free quarks, which becomes strongly interacting at 
lower temperatures T = (1 — 3) Tc, see e.g. [i^ HI] for 
recent reviews and a complete set of references. 

In Ref. [l^l the following parametrisation was proposed 
for the contribution of the coloured states to the sub- 
tracted pressure APc = PciT, n) - Pc{T, 0) (we slightly 
simplify the notation): 



APc 

y4 



(F,(r)) [cosHpLu/T) + cosh{^ld/T)) 
+Fgg{T){cosh{2nu/T) + cosh(2/Xd/r) 

+ cosh{i^iu + ^^d)/T)). (19) 



The susceptibilities at zero chemical potential can be eas- 
ily computed from Eq. (|19p . and we recognise that their 
ratios allow the identifications of the relevant degrees of 
freedom. These prediction for the susceptibilities ratio 
was contrasted with the numerical results, finding a poor 
agreement. 

The imaginary chemical potential approach gives the 
possibility to check directly the consistence of various 
phenomenological models by analytically continuing from 
real to imaginary /x. We can subject our data to the same 
analysis by analytically continuing Eq. from real to 
imaginary chemical potential. Setting ^isospin — Mm " 
/id — 0, and including the contribution from baryons and 
tetraquarks we get: 

AP 

^ = F,(r)cos(A*/T)) + F,,(r)cos(2/./T) (20) 



+ Fggg (T) C0S{3 ^i / T) + + Fgggg (T) C0S( ^/T) 



giving in turn: 

n{nj,T) = Fg{T)sm{nj/T) + 2Fgg{T)sm{2fii/T) (21) 
+ 3Fygg{T) sm{3^ii/T) + 4Fg,,(r) sin(3M//r) 

From the point of view of the imaginary chemical po- 
tential analysis, checking these forms correspond to per- 
form a Fourier analysis of our results. We have then fitted 
our data to the form 



PK{^^I) = Y!"^! sin(j//,/r) 



(22) 



The results of the fits are shown in Figure [TOl The re- 
duced ranges from 84 to 2.85 {Fi to F4) but the errors 
on the parameters grow big and the parameters them- 
selves are not stable. We summarise the results in Table 
1, and we conclude that, even if the trigonometric fits 
might eventually converge, it is hard to attach any simple 
physical interpretation of the parameters ^1,^2,^3,^4 
as contribution from free quarks, diquarks, baryons and 
tetraquarks. 

This result is not unexpected, as the Fourier parametri- 
sation Eq. (|^^ is not compatible with the critical fits 
Eq. Does this mean that the occurrence of coloured 
bound states is ruled out? 
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FIG. 10: Results of the fits of n{fj,i) to the trigonometric 
functions Eq. (|22|l . see text for details 

TABLE I: Parameters of the Trigonometric Fits 





F2 




Fi 


X^/d.o.f. 


-0.110(1) - - - 84 


-0.071(3) 


-0.023(2) 






11.11 


0.028(15) 


- 0.114(14) 


0.029(4) 




4.18 


0.43(11) 


-0.55(13) 


0.257(66) 


-0.049(14) 


2.85 



Not really: apart from the fact that extending an 
Hadron Resonance Gas description to coloured states is 
anyway a non-trivial assumption, given the non-trivial 
interactions that one could expect. Consider that the 
masses themselves should depend on the chemical poten- 
tial ,6,, .45j] - the so-called BKS effect. From the perspec- 
tive of the present study, the BKS effect is indeed very 
natural, in view of the phase diagram in Figure [TJ and 
the related analysis of the data in term of critical be- 
haviour which we have presented above. Remember, in 
fact, that the coefficients of the Hadron Gas parametri- 
sation in Eq above represent the contribution of the 
spectrum of resonances, hence their temperature depen- 
dence is inferred from the one of the masses. 

For instance, in Ref. [1] it was proposed that 



coloured 



n.5Te ((T/3T,)0-5 + Q.m/{T - Tj) (23) 



yielding the decoupling of the coloured masses at the crit- 
ical point. A similar decoupling should take place at 
the critical line at negative /i^, which is not compatible 
with the simple factorisation of the terms depending of 
temperature and fugacities implied by the HRG model 
EqdUl). 

In Ref. [45;] it was underscored that if the derivatives 
of the masses with respect to the chemical potential 



M"(T) 



d^M{T, /i) 



(r,/i = o) 



(24) 



chemical potential susceptibilities as probes of particle 
contents has to be revised, and, more generally, the de- 
coupling of the prefactor and the simple trigonometric 
factors predicted by the HRG model is no longer true. 
Hence, we cannot use simple forms as those of Eq (fT9)) 
to assess the particle content of the sQGP gas. 

Our analysis supports this point of view. It is inter- 
esting to note that indeed a recent direct calculation [5l| 
of the coloured spectrum indicates the survival of heavy 
coloured states above Tn- 



VI. SUMMARY AND OUTLOOK 

We have studied the critical behaviour of the system 
in proximity of the critical endpoint of the chiral and 
RW line in the negative /i^. We have given a simple 
description of the non-perturbative features of the sQGP 
phase, based on the analysis of the critical behaviour in 
the imaginary ^ plane. We have proposed an EOS of the 
form 



with a ~ 0.3, which accounts nicely for all the features 
of the numerical data. The exponent would read a = 1 
for a Stefann-Boltzmann-like law. 

From a more mathematical point of view, the proposed 
parametrisation is a Pade' approximant of order [2,1], as 
appropriate in the standard application to critical phe- 
nomena. The results thus obtained can be, in principle, 
analytically continued within the entire analyticity do- 
main. Practical limitations - discussed at length in pre- 
vious work - do arise because of numerical accuracy. 

As for the particle content of the system, our results 
suggest that a fit to 



AP 

ji4 



= Fq{T) cos(Ai/T)) + Fqq{T) cos{2fi/T) 

+ Fggg{T) COfi{3 / T) + +Fgggg{T) C0S{4pL/T) 



are large enough, the simple interpretation of the zero 



cannot afford any definite conclusion. This does not come 
as a surprise. The masses themselves, hence the coeffi- 
cients, will depend on fj, in some complicated way, which 
should anyhow conjure to give the simple behaviour ob- 
served in the data. The fact that the masses themselves 
should depend on the chemical potential while approach- 
ing the critical endpoint offers a simple realisation of the 
BKS mechanism. It should also pointed out that ex- 
tending a Hadron Resonance Gas description to coloured 
states is anyway a non-trivial assumption, given the non- 
trivial interactions that one could expect. 

It would be very interesting to confront the numerical 
results in a broader range of temperatures with these 
ideas, as well as with analytic calculations and phe- 
nomenological models. Future work should hopefully be 
able to give a coherent account of critical behaviour, high 
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temperature expansions and particle contents in the re- 
gion of the strongly interactive Quark Gluon Plasma. We 
hope that the simple description offered here might be of 
help in building such a complete picture. 
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